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Abstract 

The redistribution of intrinsic quantum noise in the quadratures of the field generated 
in a sub-threshold degenerate optical parametric oscillator exhibits interesting dependences 
on the individual output mirror transmittances, when they are included exactly. We present 
here a physical picture of this problem, based on mirror boundary conditions, which is 
valid for arbitrary transmittances and so applies uniformly to all values of the cavity Q 
factor representing in the opposite extremes perfect oscillator and amplifier configurations. 
Beginning with a classical second-harmonic pump, we shall generalize our analysis to apply 
to finite amplitude and phase fluctuations of the pump. 


1 Introduction 

A degenerate optical parametric oscillator (DOPO) has long been considered a nearly ideal squeez- 
ing device when operated just below threshold. The quantum fluctuations of the generated 
sub- harmonic field are rather immune to spontaneous emission since the two- photon transition 
governing the parametric down-conversion process sees no resonant intermediate levels. 

Nearly all prior work dealing with this problem [1,2,3] has been limited to the situation in 
which the DOPO cavity is nearly perfect. In a general approach [4,5] developed recently by the 
author and Abbott, which is based on the exact treatment of mirror boundary conditions, it 
has become possible to discuss cavity problems in quantum optics for the entire range of cavity 
transmissions possible. In the present DOPO context, this approach thus permits the extreme 
limits of a single-pass amplifier (cavity transmission —>100%) and of a nearly perfect DOPO cavity 
(cavity transmission — >0%), and all intermediate-Q oscillator configurations to be treated on the 
same footing. By employing this viewpoint (which may be viewed as a generalization of Collett and 
Gardiner’s approach [2]), we also develop a physically insightful picture of the general squeezing 

‘Permanent address 


35 PRECEDING PAGE 


blank not filmed 


problem, one which emphasizes the correlations of the input, output, and intracavity fields that 
govern the relationship of intracavity and output field fluctuations. Any reference to modes is 
altogether avoided here. 

After treating the DOPO problem with a perfectly monochromatic pump, we shall model real- 
istic experiments in which the pump field has finite amplitude and phase fluctuations. Although 
any amplitude noise of the pump has a relatively minor impact on the squeezing of the sub- 
harmonic signal field, pump phase diffusion even when it is tracked can cause a severe degradation 
of that squeezing. More detailed discussions of this problem will appear elsewhere [6]. 

2 Mathematical Formalism 

A description of the problem at hand that covers the whole gamut of cavity transmission factors 
is necessarily multimode in character. We avoid all reference to cavity modes by writing the 
fully quantized signal field inside the cavity in terms of its rightward (positive— z) and leftward 
(negative— z) propagating parts. For the positive-frequency part, this decomposition is written in 
the Heisenberg picture (HP) as 

E {+) (z,i) = (e+(z,t) + e. (*,!)«-“•*) «-*•', (1) 

in which the operators e±(z,t) have expectation values that are assumed slowly varying in space 
and time on the scale of the central wavelength 2r/ko and period 2ir/fl 0 - 

The parametric interaction of with an intense quasimonochromatic is described via 

the interaction Hamiltonian (also written in HP) in a cavity of length l filled with the parametric 
medium: 


H'dopo = Y X (J) c^ mp J* [e+ (z,l) + cl (z, <)] dz + Hermitian Conjugate (2) 

The complex pump amplitude epump is at most slowly varying in time. The constants A and x^ 
are the cross-sectional area of the cavity and nonlinear susceptibility, respectively. The notation 
used is the same as in Ref. [7]. We may write the equations of propagation for e± (z, t) in the 
slowly- varying envelope approximation as 

in which the nonlinear polarization waves p± L (z,t) driving the parametric interaction are given 
by a functional differentiation of the quadratic interaction Hamiltonian (2): 


p± L {*, 0 = m) 

Thus, on combining (3) and (4), we have the following generalization of the single-mode equations 
describing the parametric amplification process: 


(a* ± cdt ) ~ ±9e * (z ’^’ 


(5) 
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Fig. 1. The DOPO Cavity with End Mirrors at z = 0 and z = l. 


To complete the formalism, we supplement Eq. (5) with boundary connections of the intra- 
cavity e±(z, t) fields with the input vacuum fields. These connections are 

e+(0, t ) = -fe_(0, f) + fe+“( °> 0? 

e.{l,t) = -f'e + (*,0 + fV“(*,<), 

in which e.^ are the two traveling pieces of the vacuum field entering the cavity through its 
mirrors at z = 0 and 2 — 1 with inside- to-outside reflection and transmission coefficients (— r,f) 
and (— f',f') respectively (see Fig. 1). 

3 The Parametric Amplifier Problem 

Without the cavity mirrors, the oscillator reduces to the amplifier configuration in which the two 
traveling parts e + and e_ axe not coupled to each other. We may therefore concentrate on only 
one of them, say the e+— field. 

Furthermore, for simplicity, we shall assume in this section that the pump has no amplitude 
and phase randomness, so that it is strictly monochromatic. For this case, one may assume without 
any loss of generality that q is read and positive, for any constant nonzero phase <j> t of q may be 
scaled out by redefining e + (z, t) to carry a constant phase factor exp (»^,/2): 

e + (z, t) —* e + (z, t) ( 7 ) 

without altering the physics. 

By adding to Eq. (5) and by subtracting from it its Hermitian conjugate, one obtains the 
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following pair of uncoupled equations for the quadratures of e + : 

(I + ;l) X * M = (I + ;|) 


( 8 ) 


where X + (z,t) = | (e+(*,t) + e+(z,f)) \Y + (z,t) = Jj (e + (z,f) — e+(z,<)) are the in-phase and 
jt/ 2 out-of-phase quadratures. The solution of Eqs. (8) is straightforward in terms of the retarded 
time variable, t = t — z/c: 


X+{z,t) = * + (<M - F+(z,t) = Y, (0 ,t- z/c)c -'* ; (9) 

which represents a phase-sensitive amplification process characteristic of the parametric interac- 
tion. These solutions are entirely equivalent to the following time-evolution equations 

X + (z,t) = X + {z - d,0)e' rf ; F + (z,t) = Y + (z - ct,0)e~'«. (10) 

The linear relationships of Eqs. (9) or (10) indicate that both the expectation value and fluctu- 
ations about it of the X + — quadrature (K + — quadrature) of the signal field amplify (attenuate) by 
the same factor. This statement, valid both classically and quantum-mechanicaUy, clearly implies 
that any noise initially present in the signal is stretched along the ^-quadrature and shrinks along 
the F-quadrature, as shown in Fig. 2. It is in this way that quadrature squeezing comes about in 
a parametric amplifier. 


4 The Parametric Oscillator Problem 

Our treatment of the .parametric oscillator builds upon the simple amplifier analysis presented in 
Sec. 3 by limiting z to lie between 0 and l and adding mirrors at z = 0 and at z = £, which serve 
to connect e + and e. and the input vacuum fields via (6). As in Sec. 3, we restrict our analysis 
here to a perfectly monochromatic pump wave for which Eqs. (9) describe the interaction of the 
e + wave with the medium. Similar relations may be written down for the quadratures of the e_ 
-field (integrated backwards from z = £): 

*-(*.<) = X. (/,(-(/-»)/£) «*«-»; , ... 

K( 2 ,!)= K (/,!-(/-*)/<:) «-•«'-» k ' 

Since we are ultimately interested in calculating the quadrature squeezing of the intracavity 
field e±(z,t), we concentrate here onwards on the quantum fluctuations alone of the various 
quadratures. We first consider what the implications of the boundary connection relations (6) 
are for the fluctuations. Since (f, t) and (f', P) are all real, these relations are formally the same 
as those obeyed by any particular quadrature of e± and e± ac fields, including their X— and Y— 
quadratures separately. Furthermore, the two fields (or their quadratures) on the right-hand side 
(RHS) of each equation in (6) are uncorrelated at any t. To see this, we note, for example, that 
the e!J.® e (0, t) field entering the z = 0 mirror contributes to the e“ e (/, t’) field only after a time 
t‘ - t — 21 1 c during which the former field makes a full round trip through the cavity. Thus, 
e_(0, f) is correlated with e+ ae (0, t — 2 t/c) which is not correlated with e!J^ c (0, t), since the vacuum 
field fluctuations are essentially ^-correlated in time. In view of this lack of correlation, we may 
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Fig. 2. The Parametric Amplification Process. The X — quadrature is amplified by 
a given factor (taken to be 2 here) while the Y — quadrature is attenuated by the same 
factor. 


write for the quantum-mechanical variance of, say, the Y — quadrature of fields at the mirrors in 
terms of the power reflection and transmission coefficients (R,T) and (J^, T') (with R = f 2 , etc.) 

(AY+(0,t) 2 ) = R(AY.(0,t) 2 ) +T(AYT‘(0,t) 3 ); (12) 

(A Y.(t, t) 2 ) = R' (A Y + (£, t) 2 ) + T' (A Y™{£, t ) 2 ) , 

while setting z = £ in Eqs. (9) and z = 0 in Eqs. (11) yields for the propagation of variances 
through the medium 


(A y + (M) 2 ) = (Ay + (0,* - t/c) 2 ) (AK_(0,t) J ) = (A Y.(t,t - l/c) 2 ) e~ 2 ' e . (13) 

With the aid of Eqs. (12) and (13), we may express the retarded propagation of the Y+— variance 
at z = 0 in one complete round trip as 

(AK + (0,<) J ) = R(AY.(t,t - l/c) 2 ) e~ 2 *' + T(AY” e (0,t) 2 ) 

= t- 2 ' l R [ # ( AY^, t - tic) 2 ) + T (A Y“ e (l, t - t/c) 2 )] 

+r(Ay; BC (o,t) s ) O 4 ) 

= RR!t- A ' 1 (AV+(0, t - 2 t/c) 2 ) + t (ay?“(o, <) 2 ) 

+RT' e -w (a r. v8e (f, t - e/c) 2 ) . 
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Fig. 3. Round Trip Evolution of Fields and Their Variances. 


The foregoing sequence of mathematical steps in arriving at the round trip propagation of variances 
is shown diagrammatically in Fig. 3 to bring out the underlying physical picture. 

In steady state, the quantum statistical properties of the field do not change from one round 
trip to the next. In this long-time limit, suppressing the time entry of each variance in Eq. (14), 
we get 

[r(AK r( 0)’) + RTV«(Ay~W)] 

\ An(0) > - (1 - ft *'.-*') ’ (15) 

a result that is uniformly valid for all values of (R, T ) and (R', T') pairs (with the obvious energy- 
conservation constraints, R + T = R' + T 1 = 1). It is also worth noting that in the derivation 
of (15), the only property of the input fields used was their white-noise (^-correlated) cha racter. 
Thus, (15) applies not Just to vacuum- field inputs, Tut to arbitrary white-noise input fields. 

In the good-cavity limit, R, ft ~ 1, qt ~ 0, we recover the result of Collett and Gardiner 
generalized to allow for arbitrary white-noise input fields at the two mirrors: 


(An(0) J ) a 


t (ay;* £ (o) j ) + r (ay"“(/) j ) 

(T + T') + 4 ql 


(16) 


For vacuum-field inputs as explicitly indicated in Eq. (15), since the two input fields are statisti- 
cally identical (except for their direction of propagation), we may write more simply 


where 


(AY + (0) J ) 


(T + RT'e-w) 

(1 - RR'e-*«) N ' 


( 17 ) 


N vte s (A>7*‘(0) J ) = {AY:“(0) 2 ). 
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Note that the calculation of the variance (AX+(0) J ) of the X-quadrature of the intracavity is 
entirely analogous and is given by Eq. (17) provided q is replaced by —q everywhere. 

The degree of quadrature squeezing is the ratio (AVV(0) J ) /N^ which is generally the factor 
by which two input fields with the same quadrature variance, but not necessarily vacuum fields, get 
squeezed on altering the cavity. Detailed discussions of this quantity in both textual and graphical 
forms have been presented elsewhere, where its generalization to include arbitrary relative phase 
between the two traveling components of the monochromatic pump has also been derived (5,6). 

Having discussed the intracavity field, we now present the noise characteristics of the output 
field. Like the former field, the latter field is strongly correlated with the input fields as well. 
However, unlike the former, the output field quadratures can be easily subjected to a spectral 
analysis by choosing a sufficiently narrowband local oscillator field and integrating long enough in 
a balanced homodyne setup as was done in the original experiments [8]. We shall see that it is in 
this spectral sense that the output field exhibits a very high degree of squeezing. 

The boundary connection of the output is similar to Eqs. (6). For example, the leftward- 
traveling output field at the z = 0 mirror is a linear superposition of the transmitted part of e_(0, t ) 
and the reflected part of c" e (0, t). So any quadrature of the output field, say its Yl— quadrature, 
obeys the boundary connection formula 

*W<M) = + f>7“(o,0- (18) 

However, unlike the intracavity field, we must know the full time dependence of t), not just 

of its variance, before it can be spectrally analyzed. Equivalently, as (18) shows, we must know 
how K(0, t) evolves in time. But, that is easy to write down over a complete round trip since we 
know via Eqs. (9) and (11) how the intracavity field e± interacts with the active medium in a 
single pass through it, while Eqs. (6) tell us how the input fields e£ e leak into the cavity at the 
z = 0 and z = l mirrors. The round trip evolution of Y1(0, <) turns out to be 

Y_(0, t) = ffe- l « < y_( 0, t - 2 l/e)- tfe-WYf'i 0, t - 2 t/c) 

+ tV«*y»“(/,t-*/c), 1 

which could also have been written down directly based on physical arguments presented below. 

If K_(0, t — 21 /c) is the ^-quadrature of the cavity field just before it is incident on the z = 0 
mirror from the right then after that mirror reflection a fraction if it is reflected while a fraction 
i of the input field Y+“ c (0, t — 2 t/c) is transmitted. The two waves propagate rightward through 
the medium with their Y —quadratures attenuated by factor t~ qt . They are then reflected at the 
mirror at z = f by factor — r 7 while a fraction t‘ of the second input Y’*'“(0, t — t/c) is added to the 
circulating wave. The net field then propagates a distance l leftward through the active medium, 
with its Y -quadrature attenuated further by e -,/ as a result, to become the net field, given by the 
left-hand side of Eq. (19), a time t/c later. 

A Fourier analysis of Eq. (19) is straightforward. We shall focus only on the central (zero- 
detuning) frequency component since it has the largest noise reduction. Denoting the Fourier 
transform of a function f(t) by f(6uj), we see that for 6u = 0, Eq. (19) yields 

K(0,0) [l - rr'e- 2 ’*] = -ir , e~* t Y? e (0,0) 

+{V’'VT e (*,0), 
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while Eq. (18) yields 


VUO.O) = *K(0,0) + fV^O.O). 

By eliminating K. (0,0) between these two relations and using the energy-conservation relation 
f 2 + P = 1, one may easily show that 

(f - f'e" 2 *) V^ c (0, 0) + 0) 

n.,(0,0) = i (1-rPe-V ) ’ (20) 

whose variance is related to the spectral variance of (uncorrelated) input-field quadratures. If we 
assume that the input fields have the same spectral quadrature variance at a given frequency, such 
as is surely true for vacuum-field inputs then the spectral squeezing of the output field at zero 
detuning is by the factor 


sffi(0) s (AYUW?) / (AVT(0,0) J ) = 


(f - f'e- 2 **) 3 + Pt' 2 *- 2 '* 
(1-ff'e-W) 3 


( 21 ) 


Just as for the cavity field, the ratio ^AA^O, 0) 2 ^ / ^ AAj* e (0, 0) 2 ^ for the A-quadrature is given 
by replacing q by —q everywhere in relation (21). 

It is worth noting that just below threshold rr't 2 * 1 — » 1, the A-quadrature of the output 
field at the z = 0 mirror has infinite variance in its central frequency component, while the 
corresponding Y -quadrature spectral component has a finite variance that depends on how large 
the transmission T* of the other mirror is. In particular, for T* = 0 regardless of the value of R 
(or of T), the output V-quadrature has zero spectral variance at the line center. This is a very 
surprising result, implying as it does that even in a very low -Q but single-ended cavity the output 
field is perfectly squeezed in the spectral sense, if the parametric gain is high enough to drive the 
oscillator to its oscillation threshold. A more complete discussion of the output field, including 
the bandwidth of the squeezing spectrum, may be found in Ref. [6]. 


5 Squeezing in the Presence of Pump Noise 

In a real experiment, pump noise is inevitable. Typically, the pump field has both amplitude and 
phase noise that can be described well in classical terms alone. For example, the pump amplitude 
may have a small fluctuating piece, described in Eq. (5) via a time dependent q , 

*(«).* + ««(*), (22) 

in which 6q(t) is an Ornstein-Uhlenbeck Gaussian process with zero mean and an exponentially 
decaying two-time correlation 

(6q{t) 6q(t’)) = a 0 re -r, ‘ -r *. (23) 

The pump phase noise, on the other hand, is ultimately limited by phase diffusion which is well 
described by a classical Wiener-Levy Gaussian random process with zero mean value for the time 
derivative of the diffusing phase, 6t/>(f), and its two-time correlations: 
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(ty) = 0; 6rj>{t')) = 2DS{t - 1'). (24) 

The constants 2r and 2D are the amplitude and phase-noise contributions to the total pump 
linewidth. 

Since detailed discussions of this problem have been presented elsewhere [6], we shall restrict 
our derivations here to its relatively simple but physically revealing aspects. To begin with, we 
shall take the white-noise limit, T — * oo, for the amplitude noise. In more precise terms, this is 
the limit in which Tt/c > 1. 

Since q in Eqs. (5) and (8) is time dependent, the exponentials in Eqs. (9) and (11) have 
integrals in their exponents. For example, in Eq. (9b) one must replace 

e -7* 


for a given statistical realization of 6q. This means that the Y — quadrature variance is down by 
the factor 



6i(t-z/c)di 



c ~2f/+4ao/c 


in every single pass either leftward or rightward between the two mirrors. We used the familiar 
result that for a Gaussian random variable x, 


(e*) = c< r> eH^> 

and the fact that when Tt/c » 1, 


(25) 


(*?(*) 6q{t')) ^ 2 a 0 6(t - t '), (26) 

to obtain the preceding factor. 

A recognition of the extra factor e* a ° lc by which the Y — quadrature variance is altered when 
the pump amplitude has a fluctuating piece immediately tells us that Eqs. (15) and (17) must 
also be altered accordingly. Thus, for example, Eq. (17) takes the form 

. . (T + RT'e- 2 " 1 * 4 **") 

(AFV(O) ) = 

Since or 0 > 0, the net effect of the 6-correlated pump amplitude fluctuations is to merely reduce 
the parametric attenuation of y-quadrature fluctuations thereby leading to a smaller intracavity 
squeezing. 

Although we have not discussed the opposite, static pump amplitude noise limit, Tt/c <C. 1, 
it can be seen by physical arguments that for a given amplitude noise (6q 2 )^ , the static case 
compromises intracavity squeezing more dramatically than the white-noise case, for it is roughly 
the zero-frequency Fourier component of the pump noise spectrum that controls the steady state 
characteristics of the signal field. As the noise bandwidth T increases, a fixed amount of amplitude 
noise is partitioned into more and more Fourier components, so that the zero-frequency component 
(like any other) goes down. 
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We turn now to the computation of spectral squeezing of the output field in the presence of a 
<5- correlated pump amplitude noise. This task is quite involved when compared with the derivation 
of the preceding intracavity variance formula. One must begin with the fluctuating analog of (19) 
which may be shown to be 

Y. (0, t) = t - 2l/c)~ ff' e -"<‘)-’K«-</e)y««( 0i t - 2 t/c) 

+ (/, t — t/c) , 1 ’ 

in which 


7(0 = / [«> + *?(* - */c)\dz. ( 28 ) 

JO 

A direct Fourier transform of Eq. (27) is not possible. We must compute first the two-time 
correlation functions (K(0 t «)K.(0,O>. (K.(0,0*r*(°.O). *** (yr(0,t)y~ e (0,t')) that enter 
the output autocorrelation function (l««<(0t 0%«t(0t O) v * a Eq. (18). A Fourier transform of the 
output correlation then furnishes the spectral variance. To compute the former two correlation 
functions, we solve Eq. (27) for K.(0, t) iteratively in terms of at successively earlier times, 
one differing from the next by the roundtrip time 2 t/c: 


Y. (0, /) = -tf £ (ff') n e-’»«+’<‘>y»~(0, f - 2/(n + 1 )/c) 

n=0 

+ ? £ (rf') n t _ /(2n + l)/c) 

n=0 

in which 


(29) 


7 p (0 = r * (90 + *9 (t - z/c)] dz. (30) 

JQ 

We may use the identity (25) and the white-noise approximation (26) to obtain the useful formula 

( e -^) e -vC>)= C -(^P'W C ^ 4 ^) + ( 4 ,, J') +J ^ V> ] ( 31 ) 

— c ~(P+p')fo< e < *« e *(» , +> , ' +J » ,< ) i 

in which p < is the smaller of (p,j/). 

When combined with the 5-correlated nature of the vacuum fields, relation (29) enables one 
to secure the needed correlations from which the following output quadrature autocorrelation 
function is obtained [6]: 

(YU 0, t) YU 0, *')) = £ [P (P + 

+ £ £.(rfT +n ’^ (t-r- 2 („-n')^/c)e- 2(n+n ' ),#< e 2o « rf < n>+3n< ) ( 32 ) 

n=0 n'sO 

- p £ (ff') n IS (t-f- 2 nt/c) + S (t - f + 2n*/c)] e~ Jn «> /+Jo « nc/ + $ (f - *')] 

nxO 

in which 


(y; ae (o, o y; ac (o, o) = { y. c (<, 0 vr* (*, O) = c* (t - 0 . , 

i 

i 

I 
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As before, we are only interested here in the central frequency component of the quadrature 
spectrum. This is obtained from Eq. (32) by integrating it over (t — t') in the range (— 00 , 00 ), 
which is a trivial task due to the presence of a 6-function in every term. The resulting infinite 
sums are related to the geomteric series and can be carried out in closed form. The net result of 
these straightforward steps is the following noise reduction factor at line center: 



P (p + f 2 f /2 e~ ^ w ~ aa>c ^) (l + ff / e~ 3(t, ~°* c ^) 

(1 — fJp3 c -<(»-*«oeK) — fpe-Slw-wK) 


2P 

1 _ -«««)< 


(33) 


When the pump fluctuations are absent (ao = 0), this expression naturally reduces to result (21). 
In general, however, a graphical presentation of (33) is imperative for physical insights. This is 
done in Fig. 4 for a symmetric cavity (R = Rf). It is no surprise that as the pump amplitude 
fluctuations increase in strength, the amount of squeezing reduces for any fixed value of R (i.e. 
along a vertical line on the figure). For a fixed fluctuation strength, on the other hand, the higher 
its value the slower the squeezing increases, with increasing R, to its maximum value at oscillation 
threshold. 

A reduction of the amplitude-noise bandwidth, so that Tl/c is no longer large compared to 
1, leads to reduced output squeezing for the same reasons as for the intracavity field. It is 
worth noting that amplitude noise, being essentially multiplicative in nature (see Eq. (5)), is 
less important than pump phase noise which unavoidably couples the squeezed quadrature to the 
highly fluctuating quadrature, thereby seriously undermining squeezing. 


6 Pump Phase Fluctuations 

Even the quietest pump, such as one generated by a highly stable laser, has intrinsic random 
phase diffusion arising from the purely quantum mechanical process of spontaneous emission. 
This means that squeezing in the sub-harmonic signal field when measured relative to a fixed (or 
independently fluctuating) phase will show a time-dependent behavior as both the squeezed and 
unsqueezed orthogonal quadratures with phases slaved to the pump mix. However, if both the 
local oscillator (LO) and pump are derived from the same laser, then the reference LO phase and 
the phase of the ideally squeezed quadrature track each other. In spite of this phase tracking, 
there is a residual effect on squeezing, due to the time dependence of the pump phase diffusion 
[9], which we consider here. 

In the presence of a finite ty(i), as described by a Wiener- Levy Gaussian random process with 
moments (24), Eq. (5) has q replaced by qe^*\ and the signal quadratures X±(z,t) and Y±(z,t) 
are defined relative to the phase 6rl?{t)/ 2: 

X ± (z,t) = i [e + (r,Oe- ,f * (0/J + e-(M)e* M ‘ ,/a ] ; 

Kt(*, t) = £ [e+(z, «) e -**W a - e.(z, <)e <f *‘>' 3 ] . 

These quadratures evolve according to the matrix equation 


(34) 

(35) 
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Fig. 4. Squeezing of the Central Frequency Component of the Output Field Quadra- 
ture in a Symmetric Cavity. The full, dashed, and dotted curves represent values of 
the fluctuation parameter aoci equal to 0, 0.005, and 0.01, respectively, while the 
roundtrip gain coefficient q$l is 0.05 in each case. 


in which the column vector V±(z,t) is (JV ± (z,f), V±(z,t)) r and the <r’s are the Pauli matrices 

"'■(l o) ! *"(< •*)’ **"(• -0‘ <36) 

Although Eq. (35) is a first order equation, it is a matrix equation with the coefficient matrix 
on the RHS at any time not committing with itself at another time. This renders the solution a 
formal one in terms of time-ordered or path-ordered exponentials. The path-ordering (or time- 
ordering) has however the advantage that successive path-ordered (or time ordered) exponentials 
from one roundtrip to the next may be easily multiplied. One first combines the solution of Eq. 
(35) with the boundary connections (9) to determine the single roundtrip evolution of V+(0, t) 
to obtain a matrix analog of Eg. (19). Iterative processing of such equation leads to a formal 
solution that can, via the simplicity of writing products of time (or path) ordered exponentials 
with contiguous limits as single time (path) ordered exponentials over the entire time (or path) 
interval, be expressed in the form 

V + (0,i) = £(ff') n C(0,2/,n;t)e I * <n '*ir v “(t-2n//c). (37) 

nsO 

In Eq. (37), W vae is a column vector related to the quadratures of the two known input fields 
and C(0, 2 In; t ), a path-ordered matrix exponential involving an integral over represents the 

residua] effect of pump phase diffusion over signal noise. 

In Ref. [6], solution (37) serves as the starting point for computing the various variances and 
correlations needed for determining the steady-state intracavity quadrature variances and output- 
field quadrature noise spectrum. Eq. (37) is sufficiently complex that a statistical averaging over 
the phase noise £ri>, in spite of its Gaussian and ^-correlated nature, cannot be exactly performed in 
the involved integrals. One must settle for a series expansion of intracavity and output squeezing 
in powers of the phase diffusion constant D, which has been determined to 0(D 3 ) [6]. We refer 
the interested reader to that reference for more details. It suffices here to state that pump phase 
diffusion seems to be most important near threshold where the fluctuations in the X —quadrature 
of the cavity field have a highly slowed relaxation rate. 

7 Conclusions 

We have presented here an analysis of squeezing in a degenerate parametric oscillator that lends 
itself to an easy physical interpretation for the most part. For completeness, we have also summa- 
rized the impact of pump amplitude and phase noises of sorts encountered in a real experiment 
on the observed degrees of cavity and output squeezing. An exact analysis for the case of a finite 
pump-phase diffusion noise 64>(t) is beset by the difficulties of computing the statistical averages 
of path-ordered integrals involving Stp(t). 
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